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ABSTRACT 

Shocks around clusters of galaxies accelerate electrons which upscatter the Cos¬ 
mic Microwave Background photons to higher-energies. We use an analytical model to 
calculate this inverse Compton (IC) emission, taking into account the effects of addi¬ 
tional energy losses via synchrotron and Coulomb scattering. We find that the surface 
brightness of the optical IC emission increases with redshift and halo mass. The IC 
emission surface brightness, 32-34 mag arcsec“^, for massive clusters is potentially 
detectable by the newly developed Dragonfly Telephoto Array. 

Key words: acceleration of particles — galaxies: clusters: general — intergalactic 
medium — radiation mechanisms: nonthermal 


1 INTRODUCTION 


According to the standard model of hierarchical structure 
formation, accretion sh ocks occur around the virial radii 
of ma ssive clusters (e.g., IRvu et ^l2003l : |Schaal fc Soringell 
I2OI5I. and references therein). In these shocks, electrons 
are expected to b e accelerated by first-ord er Fermi mech¬ 
anism (see, e.g., iBrunetti fc Jon3 201^ for recent re- 


view) to have po wer-law distribution ( Blandford fc Qstrikerl 


1 19781: Ib^ 1 197811 Turbulence in the inte rgalactic medium 
(IGM; IRvu et aB l2008l : iTakizaw^ l2008l : iMiniatil l2015ll is 
another source of particle acceleration v i a second-order 


Fermi mechanism (Schlickeiser et al.l 19871 : IBrunetti et al.l 


I 2 OOII : IPetrosiajil I 2 OOII : Fuiita et al.l 1200^ 1. Relativistic elec¬ 


trons give rise to inverse Compton (IC) emission by up- 
scattering Cosmic Microwave Background (CMB) photons 
as well as radio synchrotron emission. So far, IC emis- 


ray bands both theoretically and observationally (ISarazirJ 

I 999 I: Loeb & Waxman 

2 OOOI: Totani & Kitavamal 2000l: 

Takizawa & Naitc 

200 c; 

Fuiita & Sarazin 2001: Takizawal 

2 OO 2 I: Keshet et alJ 

2003, 

l2004al. 2 OI 2 I: Petrosian et al.ll2008l: 

Kushnir & Waxman 2010l: Bartels et al. 

2015). Howeyer, 

currently there are only upper limits 

Ota et all I20l4 


ICastaldello et aP l2015l . and referenc es therein) and a few 
claimed detections of hard X-rays dReohaeli et al.l 20081) . 


but no det ection in the gamma-ray band 1 Ackermann et al] 


[ 2 ^,[ 2 ^, which may constrain particle acceleration m ech- 
anisms (e.g., IZandanel fc An^l2014l : IVazza et al.ll2015l ). Al¬ 
though IC emission from individual clusters has not yet 


* E-mail: ryo@phys.aoyama.ac.jp 


been detected, the cumulative emission from all of them 
may contribute to extragalactic ga mma-ray background 
dboeb fc Waxmanll2000l : lMiniatill2002l ). 

In this paper, we focus on the IC emission namely 
in the optical band. So far, it has been thought that 
the optical IC emission is too dim to be detectable (e.g., 
ISaraziniri999l : lFuiita fc Sarazinll2001^ . However, an advanced 
technique has recently bee n developed in the form of the 
Dragonfly Telephoto Arr ay dAbraham fc vaii nokkumll2ni4 
Ivan Dokkum et al.ll2014l ). which is optimized for the detec¬ 
tion of extended ultra low surface brightness structures and 
is capable of imaging extended structures to surface bright¬ 
ness levels below 32 mag arcsec“^ in the SDSS g band with 
a reasonable exposure time. This newly developed technique 
provides a new motivation for calculating the brightness of 
the optical IC emission in detail. Since no firm detection of 
IC emission at any wavelength has been reported as of yet, 
the optical telescopes hold the potential to bring the first 
clear detection of the IC emission from large-scale shocks 
around clusters of galaxies, which is also the first evidence 
of nonthermal processes at accretion shocks. 


We construct a simple one zone, analytical model for IC 
emission from cluster shocks, allowing us to capture the es¬ 
sential physical details as well as the parameter dependence 
of the results. Because the electrons emitting the optical 
IC emission have the Lorentz factor 50, t hey are poten¬ 
tially affected by Cou lomb energy losses fe.g.. lSarazi ni l 19991 : 
IPetrosian et al.l 20081 ). However, as we show later, this effect 
is not significant. Our model applies also to other sources 
of nonthermal electrons. We focus on IC from primarily ac¬ 
celerated electrons for simplici ty, although secondary elec- 
trons may also contribute fe.g.. lBlasi fc Colafrancescdll999l : 
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2 Yamazaki and Loeb 


iMiniatilbood : llnoue et al.ll2005| : iKushnir fc Wajcmanll2009ll . 

For the small magnetic field strength expected in clus¬ 
ters, the contribution of synchrotron emission to the optical 
brightness is negligible. 

Our paper is organized as follows. In section 2, we de¬ 
scribe our analytical model for calculating the IC emission. 
The surface brightness in the SDSS g-band is then calculated 
for fiducial parameters in section 3. Finally, we summarize 
our results and predictions for specific clusters in section 4. 
Throughout the paper, we assume a flat ACDM universe 
with cosmological parameters, h = 0.6774, = 0.3089, 

Qa = 0.6911, and fib = 0.0486 llAde et al.ll20lil 


2 ANALYTICAL MODEL OF IC SPECTRUM 

2.1 Physical quantities of cluster shocks 

A halo of mass M collapsing at redshift 2 ; has a virial radius 
Tvir, within which the mean density is Ac times the critical 
density p c(z), and circular veloc i ty at the virial radius 14 
given by (iBrvan fc Normaiill 19981 : fearkana fc Loebl[200ll ~). 


3M 


1/3 


47rAcPc(2:) ^ 

= 0.79 {I + z)~^Mpc , 


( 1 ) 


14 = 


GM 

rvir 


= 7.4 X 10^/i^/®Mi/'^w(2)^/®(l-h2)^/^cms“^ , (2) 

where M 14 = (M/10^'*Mq). The function w{z) is given by 


w{z) 


Hm Ac 
187r2 


where 

_ flm(l + Z)^ 

^ ^ + z)3 + Qa ’ 

Ac = 187r^ -h 82(111 - 1) - 39(111 


1)^ 


( 3 ) 

( 4 ) 

( 5 ) 


Note that w{z) is a monotonically decreasing function of z, 
starting at u)(0) = 0.58, decaying through ui(l) = 0.35 and 
asymptotically approaching w(z) —>■ Hm = 0.31 as 2 —>■ 00. 

For simplicity, we assume that a spherical virial shock 
is formed at rvir, and that accretion is smooth and not as¬ 
sociated with mergers of sub-units. Recent numerical sim¬ 
ulations have shown that accretion shocks are deformed 


and far from spherical (e.g., Rvu et al.l 

20031: iLau et al.l 

l2015l:ISchaal & SDringelll2015l: iNelson et al 

2 OI 5 II. however. 


emsemble-averaged gas profile shows that virial shocks exist 
near the virial radius. Future extensions of this work can be 
based on numerical simulations of non-spherical configura¬ 
tions. The ensemb le-averaged mass accretion rate onto the 
halo is written as llWhitelll994l '). 

14 ® 

M = /acc-^ , (6) 

and the shock temperature is given by, 

kT = fTpmpV^ , (7) 


where prup is the average mass of a particle (including elec¬ 
trons), and we adopt p = 0.6. The factors /acc and /t are 


dimensionless numbers of order unity. The gas density just 
in front of the shock can be written as, 

^ _ (llb/llm)M 

= 8.0 X 10“®/accw(2)(l -I- 2 )® cm“® . (8) 


This estimate i s roug hly consistent with the results given by 
IPatei fc LoebI (l2015l i if /acc « 0.5. The gas is compressed at 
the shock with shock compression ratio r, which is somewhat 
uncertain. Recent numerical simulations have shown that 
the shocks are no t so strong and their typical Mach numbe r 
is around a few (IRvu et al.l l2003l : ISchaal fc Sprint l2015l l , 
i.e., r < 4. The IGM is pre-heated and the shocks are not 
so strong at present epoch 2 ~ 0, however, at high redshifts 
2 > 1 when the IGM is cold, the shocks around clusters 
are likely to be strong (r « 4). Fortunately, we will see in 
section 3 that the IC flux does not depend on r for fiducial 
parameters. 

The magnetic field around t he shock is amplified 
as in supernova r e mnants (e.g., IVink fc Lamind I 2 OO 3 I : 
iBamba et al.l l2003l . I2005alf9l . Assuming that the en¬ 
ergy density of the downstream magnetic field consti¬ 
tutes a fraction of the downstream thermal en- 
ergy density, we estimate the magnetic field strength as 
dWaxma n fc Loe'9l2000l: i Keshet et al.ll2004'bl : iFuiita fc Katol 


I 2 OO 5 I: iKushnir fc WaxrnaJ 2 


B = (127rr^Bng/TMmp)^/^K 

= 0.72 (r4/T/acc^B,-2)^/"Mj^f w( 2 )"/®(l + zf pG(^) 


where Cb ,-2 = (^b/ 0.01) and r 4 = (?’/4). 


2.2 Injected electron spectrum 

We assume a single power-law form of injected electrons 
iV(7) = iVo7“^ , (7min < 7 < 7max) , (10) 


with a constant normalization No and a spectral index p. 

The maximum Lorentz factor, 7max, is determined by 
the balance of acceleration ti me and the cooling time. The 
acceleration time is given by llPrur 3H), 


tacc ( 7 ) 


r(r -I- 1) Tjg-fmsC^ 
r - 1 


( 11 ) 


where the gyro-factor rjg is of order unity, and the shock ve¬ 
locity that is measured in the rest frame of the shock, Ush, 
is related to 14 through Ush = [r/(r — 1)] 14. Here we as¬ 
sume Bohm diffusion with no change in the diffusion prop¬ 
erties across the shock, and no shock modification due to 
accelerated particles. We equate face to the cooling time via 
synchrotron and IG emission, 

fic/syn( 7 ) = (^2-1^2 -^ ’ (^ 2 ) 

crT(H^ -I-Hcmb)7 

wher e Bomb = 3.24 (1-1-2)^ pG, to obtain l|Loeb fc WaxmanI 


7max = 3.4 x10M/"(^)'^' 

„ (/T/accgB,-2)^^"‘Mb/^ W(2)^/^ 

?7y"[l + (B/BcMB)2]l/2 (1 +2)1/2 
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Optical IC emission from clusters 3 


based on Eqs. m and We note that the effect of syn¬ 
chrotron cooling is negligible since, 


B 


(scMB J 




= 4.9 X 10 rifrfacc^B, 




(14) 


is always small for our adopted parameters. 

In order to determine both the normalization constant 
A^o and the minimum Lorentz factor 7 min, we make two as¬ 
sumptions. One is that the production rate of the accelerated 
electrons is a fraction r\^ of the particle number input rate 
across the virial shock, Wn = (Q.hl^-m)Mj The other 
is that a fraction of thermal shock energy {3/2)kTNin 
is carried by relativistic electrons. These conditions can be 
written as 


N{'y)d'y = rj^Nir 


^7min 

/*'Ymax 


'yrUeC N{'f)d'y = Ce-feTiVu 


(15) 


(16) 


One can solve these two equations numerically for No and 
7 min, given p, 7 max, Pe and For our fiducial parameter 
set, 7 min is much smaller than 7 max and if p > 2 it is ap¬ 
proximately given by. 


7 min 


3(p - 2)kT 


= 8.4 


2 (p — l)pemeC^ 

{p — 2)Ce,-2 


(P - l)»7e.-5 




(17) 


where ie ,-2 = (^e/0.01) and ? 7 e ,_5 = (pe/lO”®). For our 
fiducial parameters, this approximate formula is accurate 
enough as long as p > 2.2. Finally, for convenience, Eq. Gi) 
can be written as. 


Nom^c^ = 0 

fip) 2 


firrip 


where 


ff J ln(7max/7min) (p = 2) 

7^7^n"[l-(7min/7max)"-^] (p / 2) 


(18) 


(19) 


2.3 Spectrum of IC emission 

Next we derive analytically the radiation spectrum of the 
IC emission for a power-law distribution of injected elec¬ 
trons given by Eq. m- A similar analysis has been done 
for s ynchrotron radiation in the stu dy of gamma-ray bursts 
(e.g., ISari. Piran. fc Naravanj[l9^ . 

The radiation power and the characteristic frequency of 
upscattered CMB photons which is ra diated from a relativis¬ 
tic el ectron with Lorentz factor 7 are (iBlumenthal fc Gouldl 

Hz]), 

P(7) = laTC7":?|^ , (20) 

/ N 4 2- 

^( 1 ) = g7 l^CMB 

= 2.05 X 10 “ 7 ^(l-b 2 ) Hz , (21) 

where hPcMB ~ 2.70 kTcMB is the mean energy o f CMB 
phot ons and we adopt Tcmb = 2.726 (1 -I- 2 ) K dFlxsenl 
I 2 OO 9 I L The spectral power, (power per unit frequency, in 


units of ergs s“^Hz“^), is proportional to v for v < v[‘y) 
and cuts off sharply at iz > yipf). The function Pv is peaked 
around 1 ^( 7 ), and its peak value is well approximated as 
~ P{‘y)/v{'y). Note that P,/,max is independent of 7 . 

The above description of P,^ is only suitable when the 
electron does not lose a significant fraction of its energy to 
radiation. This requires the characteristic cooling time of the 
electrons to be longer than the dynamical time of a cluster, 
which is given by, 

tdyn = ^ = 1.5 w{z)~^^^{l + Gyr . (22) 

yc 

Otherwise, the effect of energy loss must be considered. The 
energy loss rate of electrons is dominated by Coulomb colli¬ 
sions at l ow energies and synchrotron and IC losses at high 
energies (ISarazi rJ ll999l : IPetrosian et al.l I 2 OO 8 I I . The cooling 
time via Coulomb collisions is well approximated for rela¬ 
tivistic electrons as, 

= 3uTcIlnA '^ • (23) 

In the following, we assume that the electron density down¬ 
stream of the shock is rie « 0.5rng, where the factor 0.5 
represents the number fraction of electrons to total gas par¬ 
ticles, and the gas density Ug is given in Eq. ([S]). For sim¬ 
plicity, we fix a Coulomb logarithm at a value In A = 40. 
The synchrotron and IC cooling times have been already 
derived in Eq. m- One can find the Lorentz factors, 7 bi, 
7 b 2 and 7 b 3 , at which two of the three timescales, tic/syn( 7 ), 
tic/syn( 7 ) and tdyn, are balanced, such that. 


flC/syn(7bl) 

— ^dyn 

(24) 

tcimb(7b2) 

= flC/syn(7b2) 

(25) 

tcimb(7b3) 

— ^dyn 

(26) 


Since tic/syn oc 7 \ tcimb oc 7 and tdyn oc 7 °, 7 b 2 is always 
between 7 bi and 7 b 3 . Electrons with Lorentz factor 7 do not 
suffer significant cooling only if 7 b 3 < 7 < 7 bi ■ 

To find the spectral shape of the IC emission taking 
into account the electron cooling, we define characteristic 
frequencies as. 


^hi ~ 

y{phi) , i = l,2,3 

(27) 

and obtain. 


t'bl 

- 4.8 xlO^^ "^'1 Hz , 

[1 -b {B/BcMBYf 

(28) 

lZb2 

\ 40 / 1-b (B/Bcmb)^ 

(29) 

yhs 

= 1.8 X 10 “ r|/ic«'(2)(l + Hz . 

(30) 


These three frequencies coincide at a redshift 2 eq approxi¬ 
mately given by, 

l + 2eq«3.6(i^) , (31) 

where the term (B/Bcmb)^ is small and hence neglected 
[see Eq. (fTlll l. One can see iAjO < ^'b 2 < vsi for 2 < 2 eq 
and i/bi < Vh 2 < yho for 2 > 2eq. The spectral shapes are 
different for these two cases, and are treated separately in 
the following. Our results are summarized in Figure [T] 
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4 Yamazaki and Loeb 


2 . 3.1 1 C spectrum for Vhs < Vh2 < Vhi 

An electron with an initial Lorentz factor 7 > 7bi cools 
down to 7 bi in the dynamical time tdyn- Recall that the peak 
value of the instantaneous emissivity, Pi.,max, is independent 
of the electron energy. Thus, the average emission power at 
frequency v is proportional to the cooling time of electrons 
with Lorentz factor 7, satisfying v = *^(7) oc 7^. Therefore, 
the average spectrum scales aflP. tX flC/syn( 7 ) fX for 

r'bi < < I'ky)- For y < vm, the spectrum has a low-energy 

tail, Pv oc V. For v > v('y), the spectrum steeply decays. The 
averaged spectrum from such electrons has a peak at Vhi. 

An electron with an initial Lorentz factor 7b3 < 7 < 7bi 
does not suffer significant cooling, and the radiation spec¬ 
trum is Pv oc V for V < v{'y) with a sharp cut off for 
y > y{'y). An electron with an initial Lorentz factor 7 < 7b3 
suffer energy loss via Coulomb loss, so that we obtain 
Pi, oc tcimb(7) oc y^^^ for v < 1/(7) with sharp cutoff for 

y > y{'y)- 

To calculate the net spectrum from a power-law distri¬ 
bution of electrons, one needs to integrate over 7. There are 
three different cases, depending on 7min, in which (A) 7b3 < 
7bl 7 7min, (B) 7 b 3 < O^min < 7bi and (C) 7 min < 7 b 3 < 7 bl- 
Below we provide additional details on these regimes. 

• (A) 7 b 3 < 7 bl < 7min (i.C., t^bS < Vhl < i^m): 


In this case, all the electrons cool down to 7bi. The spec¬ 
tral power is given by. 


(A1 

yh^, = { (A2 

Lo{yjn/yhi)^^'^{v/yinf (A3 


where. 

I'm = J^(7min) . 


Vm < y) 

vhi <y < ym) , (32) 
y < Vhi) 


( 33 ) 


In this case, only those electrons with 7 > 7bi can cool. 

We have, 

r Lo(i'/i^m)^""'’F 2 (B1 : i/bi < y) 

yLi, = U{ymlymfl^(y (B2 : Vm<v< vm) , (37) 

y Lo(l'm/^'bl)^^^(i^/j^m)^ (B3 : y < ym) 

where Lo is given by Eq. 

• (C) 7min < 7b3 < 7bl (1-6., J^m < I'bS < l^bl): 


In this case, electrons with 7 > 7bi and 7 < 7b3 can 
cool. The electron spectrum has a break at 7b3, below 
which it is well approxima ted with the s tationary solution, 
oc iV(7)tcimb(7) oc 7^”^ dSarazinl Il 999 ll . Hence, we have 
Li, oc yC'-v)l‘C2 y„i < y < r'b3. Therefore, the luminosity 
is given by, 

(Cl : yhi < y) 

(C 2 : yhz <y < t'bi) 
(C 3 : y^n < y < yhz) 
(C 4 : y < Via) 

where Lo is given by Eq. 


I Lo(ymfybi)^^^(y/ym)^^~^^'^^ 

I Lo (ym/ ybl ) ^^^ (t'm/r'bs)^^^ (y/ 

I Lo(l'm/t'bl)^^^(t'm/ri,3)^^^(z^/t'm)®^^ 


2 . 3.2 IC spectrum for z/bi < r'b2 < t'bs 

In this case, the electrons suffer significant cooling through¬ 
out. Hence, an electron with an initial Lorentz factor 7 > 7b2 
cools down through 7b2 until which it produces Pi, oc 
and further lose its energy to form Pi, oc y^^^ below r'b2. 
The average spectrum from such electrons has a peak at yb2. 
Similarly, an electron with an initial Lorentz factor 7 < 7b2 
produces the average spectrum P^ oc y^^^ for y < y{'y)- 
calculating the net spectrum for power-law electron distri¬ 
bution, there are two different cases in which (D) 7b2 < 7m 
and (E) 7m < 7b2. Below we describe them in detail. 


To determine the normalization constant, Lq, one can use 
the fact that electrons with 7 > 7bi lose almost all their 
energy via IC emission, that is, the luminosity in the 
regime ym < v can be written as dLoe b Sz W axmarJ uOOOl : 


iKeshet et ahllioo^ : iKushnir fc WaxmarJ 2009 ll . 


L„dy = 'ymeC^ N ('y)d'y 


( 34 ) 


Using Eqs. (nni and GHl) and dv/y = 2 d'yl'y for iz oc 7^, we 
obtain. 


Lo = 


where. 


g{p) = 


3gefcT {Clb/Clai)M 

4g{p) fimp 


2.0 X 10^3 -^T/accg^e, ^ W {zf ^ ^ {1 + S^fe) 


ff(p) 


ln(7max/7min) (p = 2) 

^ [1 - (7min/7max)^"^] {p + 2) 


(36) 


• (B) 7 b 3 < 7 min < 7 bl ( 1 - 6 ., i^bS < i^m < t/bl): 


• (D) 7b2 < 'ymin (i'6., /^b2 ^ ^m)- 


Electrons with 7 > 7min cools to make a spectrum Liy oc 
J/-P/2 £qj- while electrons with 7min forms the spec¬ 
trum below iZm. The luminosity is given by, 

r Lo(iz/iZm)‘=-^^/= (Dl: ym<y) 

yLi, = < Lo{ylvmY^'^ (D2 : yb2 < y < ym) , (39) 

i Lo(iZm/r'b 2 )(i"/t'm)®^^ (D3 : y < yb 2 ) 

where Lo is given by Eq. dSH). 

• (E) 7 min < 7b2 (i.e., < yb 2 )'. 


In this case, electron distribution for 7min < 7 < 702 is 
proportional to iV(7)tcimb(7) oc 7^^”^. Hence, we have. 


Lo(iz/r 


x(2-p)/2 


yL, = { U(ymlyb2)(ylymi^-^'^l^ 

Lo (l'm/t'b2 ) (I'/t'm)®^^ 

where Lo is given by Eq. (| 3 ^. 


(El 

(E 2 

(E 3 


yh 2 < y) 

ym <y < yb 2 ) , (40) 
y < ym) 


^ This spectral slope is identica l to the similar case considered 
in ISari. Piran. &: Narav^ lll998h since both of the characteristic 
frequencies of the synchrotron and IC emissions are proportional 
to the square of the electron energy. 


2.4 Observed surface brightness of IC emission 

The observed surface brightness Sv (in units of 
erg s“^cm“^Hz“'^str“^) of IC emission from a cluster 
at redshift z is given by. 
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Optical IC emission from clusters 5 



Figure 1 . IC Spectrum from power-law distribution of relativistic electrons. Regimes (A) 1^3 < r'ui < t'm, (B) r'us < Vm < t^bl and 
(C) Vm < Vh3 < !A>1 are for Ub3 < r'b2 < lAl, while regimes (D) Ub2 < and (E) r-m < 1^2 are for u^i < r'b2 < iA 3 - 


(1 + z)Lv^ 

AiidL{zY^li, 


(41) 





fw icoolOO 
3 f dyn 


(42) 


where = (1 -|- z)!/, (IlIz) and 12^ are the observing fre¬ 
quency translated into the cluster rest frame, the luminosity 
distance to the cluster and the solid angle of the extended 
emission on the sky, respectively. 

The effective value of can be estimated based on 
the observed brightness profile in the sky, which depends 
on the Lorentz factor of the electrons, 'y^, — (yj^SSci x 
lO^^Hz)'^^^, emitting IC photons with observed frequency 
V [see Eq. (I21II ]. If the total cooling time of such electrons, 
tcoaxffiv) = [tic/syn( 7 v)“^ + t cimb (iv ) “^]“^ E much Smaller 
than the dynamical time of the cluster tdyn, then the emis¬ 
sion mainly originates from the region around the shock. 
The width of the emission region (in the radial direction) is 
given by the product of downstream flow velocity and the 
cooling time, <5 ~ (l/3)Ktcooi(7!/) < rvu, where a factor 1/3 
applies to the strong shock limit. The surface brightness pro¬ 
file has a rim-brightened shape. Because of projection, the 
observed apparent scale width, IT, in the sky differs from 
5. The ratio Wjd depends on the uncertain radial prohle 
of the electron dist ribution downstream of the shock, but is 
typically around 3 (iBamba et ahll^OSal l. Writing W = /w5, 
we obtain 17^ « 27r(/w/3)rvirV'ctcooi(7i/)/dAi where is the 
angular diameter distance to the cluster. On the other hand, 
if tcooi(7!^) tdyn, the cluster interior is filled with electrons 
with 7 „. Taking into account the projection effect, the sur¬ 
face brightness profile is center-Hlled. In this case, assuming 
the bright emission in the sky originates from inside the ra¬ 
dius rvir/2, it occupies a solid angle Q,i, « (7r/4)(rvir/dA)^. 
Connecting both limits, we get 


3 RESULTS 


Based on the derivations in the previous sections, we can 
now calculate the observed surface brightness Sv in the 
SDSS g-band {u = 6.3 x lO^^’lIz) as a function of redshift z 
and a halo mass M. Other parameters are hxed at the fidu¬ 
cial values, p = 2.5, r^, = Cb ,-2 = Pg = W.-s = 1, 0,-2 = 5, 
/acc = /t = 0.5 and /w = 3. The parameter is highly un¬ 


certain and depends on upstream physical quantities such as 
magnetic field and gas temperature (e.g., iMatsukivo et al] 
I2OIII: IGuo et al.l l2014all^ It can range from ~ 10 ^ to 
~ 10“^^ I Kang et al.ll2012^ . If pe is so large, then ymin is 
less than unity [see Eq. (llTIl l , so that our assumption of sin¬ 
gle power-law spectrum breaks down. Dependence of p^ is 
easily found as seen in the following. The spectral index p 
is also uncertain, however recent study of pa rticle acceler¬ 
ation at cluster shocks suggest p ss 2.25-2.5 (iKang fc Rvul 
l201ll . I2OI3I : iHong et '^12014 IGuo et ^l2014ah . In the fol- 
lowing, we also consider the cases of p = 2.0, 2.3 and 
3.0, as well as the hducial case of p = 2.5. The index p 
is related to the shock compression ratio r. For hrst-order 
Fermi acceleration, p — (r + 2)/(r — 1) in t he test-particle 
limit telandford fc OstrikerllToTSl : lBelllll978l l. In this case, r 
ranges between 2.5 and 4 for 2 < p < 3. Nevertheless, we 
hx r 4 = (r/4) = 1, which corresponds to the strong shock 
limit, because observed surface brightness does not depend 
on r for the parameters of interest. All the equations in this 
section do not depend on r. 

Figure [ 2 ] shows the surface brightness as a function of 
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redshift z, for a fixed halo mass Mu = (M/IO'^'^Mq) = 3. 
The red fine describes the fiducial value of p = 2.5, while 
the others are for different values of p with the other param¬ 
eters fixed at their fiducial values. Interestingly, the surface 
brightness increases with 2 . At low redshifts where the ef¬ 
fects of cosmological expansion are negligible, the brightness 
is almost constant, « 36 mag arcsec“^. As 2 increases, the 
surface brightness becomes larger by « 3 mag until 2 « 2. 
For this redshift range, one can see from Figure [21 that 
I'bs < Vm < Vs = (1-1-2);^ < Vhi, SO that the spectrum 
is in the regime B2 [see Eq. (1371) ]. Thus, we find, 

X + . (43) 

Since electron cooling is not significant, the solid angle of 
the emission is given by 0,^, « (7r/4)(rvir/dA)^, so that, 

Sv oc -I- 2 )“®L „3 

oc 

X . (44) 

Since w{z) is only weakly dependent on 2 , the surface bright¬ 
ness Sv increases with 2 following the scaling (1 -1- 2 )^''"^. 

The regime B2 ends when Va becomes larger than Vbi- 
This crossing occurs at 2 ~ 2 for our fiducial parameters. 
Thereafter, the spectrum is in regime Bl. When 2 further 
increases, the spectrum enters into regime Al, subsequently 
followed by the regime Dl. In these regimes, the luminosity 
is given by the same form [see Eqs. (|2^. (|271) and (12^1)]. 

X + , (45) 

and the solid angle is still given by ~ (7r/4)(rvir/dA)^, so 
that we have, 

X . (46) 

For our fiducial parameters, Eq. (M describes the scaling 
if 2 < 2 < 5.1. The flux increases by a factor of about 
3 in this redshift range. Note that the emergence of these 
regimes originates from rapid decreasing of r'bi oc (1 -I- 2 )“^^ 
for 2 > 1. As a result, the blue-shifted observing frequency 
Va = (1 -I- 2 );^ becomes larger than any other characteristic 
frequencies, r'm and r'bi (* = 1, 2, 3), around 2 eq. As 2 further 
increases, Vm oc ( 1 - 1 - 2 )^ together with r-bs oc ( 1 - 1 - 2 )'^ become 
large, finally exceeding Va, so that the spectrum enters the 
regime D2. There the electron cooling is so significant that 
the observed brightness profile is rim-brightened shape, so 
that Q.I, is small and Sv shows rapid increase for 2 > 6. How¬ 
ever, the abun dance of clusters at these high redshifts is e x- 
tremely small llBarkana fc Loe9l200ll : fWatson et ahilioi^ l. 

Although the parameter dependence is somewhat com¬ 
plicated, one can see that overall behavior is not so different 
from the fiducial parameter set; the brightness varies typi¬ 
cally by up to « 2-3 mag if one of parameters is changed 
with others fixed. Figure d] shows lines for the cases of 
p — 2.0, 2.3 and 3.0 with other parameters fixed as fiducial. 
The larger the p, the brighter the surface emission. The de- 





Figure 2. Observed SDSS g-band surface brightness of IC emis¬ 
sion from a cluster with M 14 = 3 (in units of mag arcsec“^) as 
a function of cluster redshift 2 . The red line is for fiducial pa¬ 
rameters (p = 2.5, r 4 = ^B,-2 = Pg = Pe,-5 = 1, Ce,-2 = 5, 
fa.cc = /t = 0-5 and /w = 3), while the green, purple and blue 
lines are for p = 2.0, 2.3 and 3.0, respectively. 

pendence on other parameters can be found from Eqs. (I44t 
and (l46ll . 

FigureUdepicts the surface brightness as a function of a 
halo mass M for fixed cluster redshifts assuming the fiducial 
parameter values. As seen in Fig. (2] the redshift dependence 
is small when 2 < 0.1. There, the surface brightness peaks 
at Mi 4 ~ 30, below which the spectrum is in regime B2 for 
0.55 < Mi 4 < 30 and in regime C2 for M 14 < 0.55. Func¬ 
tional forms of the brightness for B2 and C2 are given by 
Eq. (144 II . with Si, oc . Eor M 14 > 30 the spectrum 

is in regime B3, at which we have. 

Si, oc faccfT^^7^VeMr4W{z)v , (47) 

so that the brightness decreases with M. Eor 2 = 0.5, 
the peak shifts to lower masses, and larger values. Note 
that Eq. dig implies that the brightness in regime B3 very 
weakly depends on 2. The behavior for 2 = 1.3 is the same 
as for lower redshifts as long as M 14 < 33, above which, 
however, the spectrum is in the regime A3 and the surface 
brightness again increases with M and 2 according to the 
scaling, 

S„ oc /acc77eM^/®w(2)^/®(l -|- zf V . (48) 

For higher redshifts 2 = 2.0 (2 = 4.0), the spectrum enters 
regimes C2, B2 and Bl (El and Dl) in turn towards higher 
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Figure 3. Characteristic frequencies, as measured in the rest 
frame of a cluster, for our fiducial parameter values with a halo 
mass Mi 4 = 3, as a function of redshift z. The green, blue, purple 
and light blue lines indicates i^bi? ^3 and z/m, respectively. 
The red line represents Us = (1 + z)^, which is the observed 
frequency blue-shifted to the rest frame of the cluster. 

M. The brightness breaks at M 14 = 8.5 and 3.9 for z = 2.0 
and 4.0, respectively. After the break, the brightness still 
increases with M according to the scaling, 

+ . (49) 

In summary, the surface brightness follows simple scaling 
behavior, oc for the typical range of expected 

parameters. 


4 DISCUSSION 

Using a simple analytical model, we have calculated IC emis¬ 
sion in the SDSS g-band from relativistic electrons acceler¬ 
ated in galaxy clusters, taking into account the effects of 
Coulomb, synchrotron, and IC energy loss of the emitting 
electrons. For our fiducial parameters, at z < 2 and M < 
10^®Mq, the spectrum is in the regime B2 or C2, in which 
Sv oc {I + , where p is the power-law 

index of electron distribution (see Fig. 01). If the value of p 
is inferred from radio synchrotron emission, one can predict 
the spectral index of the optical IC emission. In this paper, 
we have not taken into account the possibility of reaccel¬ 
eration in the downstream turbule nce, which enhances the 
abundance of relativistic electrons llSchlickeiser et al]|l987l : 


Figure 4. Observed SDSS g-band surface brightness (in units 
of mag arcsec“^) as a function of a halo mass M for our fiducial 
parameters (p = 2.5, r 4 = ^b,-2 = Vg = he,-5 = 1, ?e,-2 = 5, 
/acc = /t = 0.5 and /w = 3). Different lines correspond to 
different redshifts z of the cluster. 


iBrunetti et al.l[200ll : |Petrosian|[200ll : lFuiita et al.ll2003ll . re¬ 
sulting in brighter IC emission. We also assume a spherical, 
virialized shock as a result of smooth accretion. If the cluster 
is not dynamically relaxed and is still being formed, e.g. as 
a result of a merger of two sub-components, then the shock 
structure would be different and turbulent energy would be 
enhanced. In this case, we expect brighter IC emission due 
to more efficient reacceleration. Hence, our estimate for the 
IC flux and surface brightness should be regarded as a con¬ 
servative lower limit. 

It would be natural to assume that electron accelera¬ 
tion occurs at the virial shocks. In merging clusters, radio 
relics have been detected as possible signatures of the merger 
shocks, although clear evidence for it had not been identi¬ 
fied as of yet. Observations at other wavelengths would be 
helpful for shock identification. Proton acceleration at the 
virial shocks could also occur, altho ugh we have only a few 
observational implications of it fe.g.. lFuiita et af]l2013h . Ob¬ 
servations of optical IC emission may have the advantage of 
enabling the identification of shocks, because the angular 
resolution of the optical telescopes is generally much bet¬ 
ter than hard X-ray or gamma-ray telescopes. According to 
the theory of diffusive shock acceleration, electrons respon¬ 
sible for the optical IC emission cannot penetrate far up¬ 
stream relative to the shock front. Hence, the IC emissivity 
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is expected to have a steep rise across the shock front. The 
detection of such a sharp jump in the IC emission would 
flag the position of the shock front, and provide evidence 
for particle acceleration there. In our spherical model, how¬ 
ever, IC-emitting electrons are not rapidly cooling in most 
cases, resulting in center-filled shape of the brightness pro¬ 
file on the sky. Thus, the profile may not have a sharp rise in 
projection. Possible exceptions might be expected for more 
realistic, non-spherical cases, where sharp feature could be 
found on the sky at the location of shocks. Detailed mor¬ 
phological studies based on numerical simulations, like the 
case of radio relics (e.g., ISkillman et al.l [20131 : iHong et al.l 
l2015l) . are needed in this case and go beyond the scope of 
this paper. 

Our model predicts the surface brightness of IC emis¬ 
sion from several massive clusters, as summarized in Ta¬ 
ble 1. Larger fluxes are expected for larger M and 2 if the 
spectrum is in regime B2 or C2 [see Eq. O]. The pre¬ 
dicted brightness ranges between « 32 and 35 mag arcsec“^. 
Such ultralow surface brightnesses should be detectable 
with the recently developed D ragonfly Telephoto Array 
ll Abraham fc van Dokkunj 120141 '). More detailed prospects 
for the observation of individual clusters will be given else¬ 
where. 

Our model also predicts the color of the IC emission, 
which is bluer than starlight. If we assume a power-law form 
of the IC emission, as Sv oc i'~°‘, then the color g — r is 
calculated as 

5 

9-^= -alog^f^{ug/vr) = Q.2<da , (50) 

where Vg = 6.3 x lO^^Hz and Ur = 4.8 x lO^^Hz are SDSS 
g-band and r-band frequencies, respectively. In the case of 
spectral regime B2 or C2 (i.e., a = (p — l)/2), we obtain 
g — r — 0.19 for p = 2.3. If the spectral regime is in B3 (i.e., 
a = —1), then g — r = —0.29. These values are distinguish¬ 
able from stellar components, such as diffuse faint emission 
of brightest cluster galaxies (BCG), satellite galaxies, and 
intra cluster light (ICL)^ whic h typically shows g — r > 0.7 
fe.g.. [Montes fc Traiillollioill . 

Gamma-ray and hard X-ray observations have given 
upper limits, which constrain our model parameters. Fig¬ 
ures [5] and [ 6 ] show that for our fiducial parameters (namely 
p = 2.5), both IC gamma-ray and X-ray fluxes for nearby 
{z <C 0.1), massive (M 14 = 30) clusters exceed current obser¬ 
vational upper limits, which are A''(0.2 — 100 GeV) < 10~^- 
10~^ photons cm“^s“^ for gamma-rays dAckermann et al] 
l2010lj a nd ^(12 - 60 k eV) < 10““-10“^°erg cm-^g-i for 
X-rays (lOta et al.ll2014h . Although the case with p < 2.5 is 
unlikely for such nearby massive clusters, our model predic¬ 
tions can be lower than these current observational upper 
limits for less massive (M14 < 10 ), higher redshift (2 > 0 . 1 ) 
or steeper electron index (p > 2.8jj. Note that we conserva¬ 
tively adopt p = 2.5 in this paper, but that higher values of 
p as implied by the gamma-ray and X-ray limits for low red- 


^ Some X-ray observations for specific clusters have pro¬ 
vided tight upper limits, such as F( 12 — 60 keV) < 5 X 
10“^^erg cm“^s“^ dNishino et ahll^OlOll for Perseus {z = 0.018 
and M 200 = 1-2 x IO^^Mq), while our model overpredicts X-ray 
flux by a factor of 3 larger even for p = 3.0. In such cases, smaller 
value of (,$ 0 . 01 ) may be required. 





Figure 5. Observed IC gamma-ray flux N{0.2 — 100 GeV), 
in units of photons cm“^ s~^, as a function of redshift 2 : for 
most massive (M 14 = 30) clusters. The red line is for fldu- 
cial parameters (p = 2.5, r 4 = ^b ,—2 = Vg — ^e ,—5 = 1, 
^e ,-2 = 5, face = It = 0.5 and /w = 3), while the purple, 
light blue and blue lines are for p = 2.3, 2.8 and 3.0, respectively. 
Mass dependence of the gamma-ray flux is given by Eq. list , 
as 7V(0.2 — 100 GeV) oc because observed frequency 

blue-shifted to the rest frame, Us, is always larger than any other 
characteristic frequencies, z/bi? ^b 3 ^m- Observed up¬ 

per limits for specific clusters lie 7V(0.2 — 100 GeV) < 10“^- 
10“®photons cm“^s“^ llAckermann et al]|2010h . 


shift clusters would make our predicted optical IC flux some¬ 
what higher (see Fig. 0. Another independent way to have 
lower gamma-ray flux is to adopt unusually large pg S> 10^. 
In this case, 7max is small enough to give spectral cutoff of 
IG emission at energy below the Fermi band. 

The IC emission from the virial shock of clusters could 
be enhanced due to relativistic electrons produced in su¬ 
pernovae that gradually diffuse out of the cluster. The lat¬ 
ter contribution depends on the star formation history and 
the diffusion time of relativistic electrons within the clus¬ 
ter. The diffusion time of electrons depends on the unknown 
configuration of magnetic fields. If the fields are radially 
aligned in the outer envelope of clusters (as expected from 
radial infall or the magneto-thermal instability; see e.g., 
[Parrish et al.l l2012h . then the diffusion time there would 
be of order the light crossing time of the outer parts of 
the cluster, i.e. millions of years. The Fermi satellite has 
placed tight limits on this contribution in cluster cores based 
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Table 1. Predicted surface brightness in SDSS g-band for specific clusters. 


Name 

z 

[1014m” 

Surface brightness^ 

[mag arcsec”^] 

p = 2.0 p = 2.5 p = 3.0 

Reference 

IDCS J1426.5-I-3508 

1.75 

5.3 

34.1 (C2) 

32.1 (B2) 

31.2 (B3) 

Stanford et al. (2012) 

SPT-CL J2106-5844 

1.132 

9.8 

34.2 (B2) 

32.1 (B2) 

31.8 (B3) 

Williamson et al. (2011^ 

ACT-CL J0102-4915 

0.870 

22.3 

33.7 (B2) 

32.1 (B3) 

32.7 (B3) 

Menanteau et al. ^20121 

SPT-CL J2344-4243 

0.596 

25.0 

34.0 (B2) 

32.3 (B3) 

32.8 (B3) 

McDonald et al. (2012) 

MS 1054-0321 

0.83 

12 

34.4 (B2) 

32.3 (B2) 

32.0 (B3) 

Jee et al. (2005) 

SPT-CL J0658-5556 

0.296 

31.2 

34.3 (B2) 

32.5 (B3) 

33.1 (B3) 

Williamson et al. (2011) 

XDCP J0044.0-2033 

1.579 

4.4 

34.5 (C2) 

32.6 (B2) 

31.7 (B2) 

Tozzi et al. (2015) 

SPT-CL J2337-5942 

0.775 

10.5 

34.7 (B2) 

32.6 (B2) 

31.9 (B3) 

Williamson et al. 120111 

Coma 

0.0232 

27.8 

35.0 (B2) 

32.7 (B2) 

32.9 (B3) 

Kubo et al. (2007) 

MACS J1206.2-0847 

0.44 

14.1 

34.9 (B2) 

32.8 (B2) 

32.2 (B3) 

Presotto et al. (2014) 

Abell 2390 

0.228 

18 

35.1 (B2) 

32.9 (B2) 

32.5 (B3) 

Carlbere et al. (1996) 

XLSSU J021744.1-034536 

1.91 


35.0 (C2) 

33.3 (B2) 

32.6 (B2) 

Mantz et al. (2014) 

Abell 2744 

0.3064 

70 

33.4 (B2) 

33.3 (B3) 

33.9 (B3) 

Montes & Truiillo 120141 


“ Values of M 200 are taken from references. 

^ For fiducial parameters ot her than p. Spectra l regime is also shown in parentheses. 

^ Since only M 500 is given in iMantz et al.l ll2014l i , we estimate M 200 assuming isothermal density distribution 
as M 200 = I. 58 M 500 . 


on the lack of Kamma- r ay em ission at 0.2-100 GeV there 
l|Ackermann et al]l2010l . |2014l . see also Vazza et al. 2015). 

Diffuse optical emission from Thomson scattering of 
starlight could be comparable to the IC emission. We 
roughly estimate the flux of the scattered light emission as, 






L.tt 

47rd_L(jt)2 


(51) 


where L* is the bolometric stellar luminosity and rx is the 
optical depth for the Thomson scattering. The typical mass- 
to-light ratio within the virial radius ry j r is g i ven by M/L ~ 
2.5 X 10^ Mq/Lq fe.g.. ISheldon et ^l200Sj : iHoIland et al.l 
l2015h . so that L« « 1.6x 10 ^®Mi 4 erg s“^. The optical depth 
tt can be estimated as, 


Tt ^e^T^vir 

- 3.3xlO-*r4U,,Mli^w{zf^^{l + zf . (52) 

Note that this estimate may be upper limit since we im¬ 
plicitly assume that most of the energy of the starlight is 
contained in the SDSS g-band in the observer frame. In Fig¬ 
ure [71 we show as a function of redshift z (green lines), 

comparing with the IC flux (red lines). One can see that the 
IC emission dominates if M 14 > 3. Since the spectrum of 
the scattered light is very different from that of the IC emis¬ 
sion, color measurements can be used to distinguish between 
them. 

The BCG and ICL could also constitute a diffuse 
background. The ICL emission has be en calculated based 
on the cosmologica l simulation (e. g., iRudick et al.l I2OIII : 
iLaporte et al.| ioj^Cui et_al.l l2014l '). According to the re- 
cent result bv icui et al. ( 2014h . the ICL brightness could be 
around 30 mag arcsec”^ at Xvir for the most massive clusters. 
On the other hand, many observations measure the bright¬ 
ness profile of the BCG- bICL only in the interior of clusters 
due t o its faintness fe.g.. |Presotto et al.ll201 j : lDeMaio et al.l 
I2OI5II . Extrapolating linearly the observed r adial profile of 
BCG-blCL of a cluster MACS J1206.2—0847 ijPresotto et al.l 


I 2 OI 4 II to its virial radius of Xvir 2.3 Mpc yields emission 
that is dimmer than ~ 35 mag arcsec”^, so that IC emission 
is brighter for this cluster (see Table 1). Even if BCG+ICL 
has comparable brightness to the IC emission, their color 
difference can be used to separate them from each other. 

IC emission also lies in infrared bands; however, the 
flux is rather small because typically Us < i^m (i.e., it is in 
regime B3). Furthermore, mid and far-infrared bands may 
be dominated by dust emission (lYamada fc Kitavamallioosi : 
iKitavama et al.fl2009l l. 

Finally, we remark on IC emission as a possible back¬ 
ground emission for other purposes. For example, as already 
discussed above, the BCG-I-ICL brightness around the virial 
radius rvir could be comparable to th e IC emission. The sur - 
face brightness of dwarf galaxies (e.g., [Herrmann et al.]l2013f l 
is also in some cases less than 32 mag arcsec”^ at the pe¬ 
riphery of the galaxies. The IC emission could confuse or 
disguise these different emissions from member galaxies in 
massive clusters. 
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Figure 6. Observed IC X-ray flux _F(12 — 60 keV), in units 
of erg cm~^ as a function of redshift z for most massive 

(Mi 4 = 30) clusters. The red line is for fiducial parameters (p = 
2.5, r 4 = |b,-2 = rig = ??e,-5 = 1, ?e,-2 = 5, /acc = /t = 0.5 
and /yy = 3), while the purple, light blue and blue lines are for 
p = 2.3, 2.8 and 3.0, respectively. Mass dependence of the gamma- 
ray flux is given by Eq. 1145 ll . as F(12 — 60 keV) oc 
because observed frequency blue-shifted to the rest frame, Us, 
is always larger than any other characteristic frequencies, Uhl, 
J^b 2 i r'hs and Um- Observed upper limits for specific clusters lie 
F(12 — 60 keV) < 10“^^-10“^'^erg cm“^s“l llOta et al.ll2014h . 
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